Chiral low-energy constants Lio and Cst from hadronic r decays 
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Using recent precise liadronic r-decay data on tlie V — A spectral function and general properties of QCD such 
as analyticity, the operator product expansion (OPE) and chiral perturbation theory (xPT), we get accurate 
values for the QCD chiral order parameters L\q and C^-j. At order we obtain L\q{Mp) = —(5.22 ±0.06) • 10~^, 
whereas at order / we get Llo{Mp) = -(4.06 ± 0.39) • 10"^ and Ch{Mp) = (4.89 ± 0.19) • 10"^ GeV"^ 



1. Introduction 

Hadronic r-decay data are a very important 
source of information, both on perturbative and 
non-perturbative QCD. Of special interest in or- 
der to study non-perturbative QCD quantities is 
the difference of the vector and axial- vector spec- 
tral functions, because in the chiral limit the cor- 
responding V—A correlator is exactly zero in per- 
turbation theory. 

The T data can be used to determine the pa- 
rameters of xPT [1]) the effective field theory of 
QCD at very low energies (a Taylor expansion in 
external momenta and quark masses). At lowest 
order, ©(p^), the SU(3) xPT Lagrangian has only 
two parameters, the pion decay constant /tt and 
the light quark condensate. At 0{p^) twelve more 
low-energy constants (LECs) appear (ii=i,...,io 
and -ffi,2), whereas at 0{p^) we have 90 (23) ad- 
ditional parameters Ci in the even (odd) intrinsic 
parity sector [2] . These LECs are related to order 
parameters of the spontaneous chiral symmetry 
breaking of QCD, and have to be determined phe- 
nomenologically or using non-perturbative tech- 
niques. For the Li couplings this has been done 
to an acceptable accuracy, but the Ci LECs are 
less well known. 

There has been a lot of recent activity to de- 
termine these chiral LECs from theory, using as 
much as possible QCD information [3-9]. This 
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strong effort is motivated by the precision re- 
quired in present phenomenological applications, 
which makes necessary to include corrections of 
0{p^). The huge number of unknown couplings 
is the major source of theoretical uncertainty. 

We present here an accurate determination of 
the xPT couphngs Liq and Cgy [10], using the 
most recent hadronic r-decay data [11]. Esti- 
mates of LiQ from r-data have been done pre- 
viously [12-14], altough our analysis is the first 
that includes the known two-loop xPT contribu- 
tions and then also the first that provides Cs?. 

2. Theoretical Framework 

The basic objects of the theoretical analysis are 
the two-point correlation functions of the non- 
strange vector and axial-vector quark currents 

IV>:lJ{q)^^jd'xe^'^^{Q\T{Jt^{x)Jr,{Q)^)\Q) 

= (-^^V+<zV)n«^(9') + 'zVni°)^(g2), (i) 

where denotes V^J^ = u^^^d and 117^756?. 
In particular we are interested in the difference 
n(s) = n.^'^^^y-' — n^'^^^\ and we will work in the 

isospin limit (m,j = to^) where ^^ulvil^) ~ ^■ 

From the analytic structure of the correlator 
n(s) in the complex s-plane and its OPE one can 
get the following two sum rules (see ref. [10] for a 
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careful derivation) 
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that represent the starting point of our work. The 
interest of these two relations stems from the fact 
that the effective parameters L^q and Cgy can 
be extracted from the data and the r.h.s can be 
rigorously calculated within xPT in terms of the 
LECs that we want to determine. From the re- 
sults of ref. [15] we get 

2fjml + n(0) = -8Ll,ifi) + GtM 

+ 16 GM + + HUp) + 0{p') , (4) 

where the functions G™^ (^) , H^^j^ (fi) are correc- 
tions of order p"^ generated at the n-loops level. 
We omit their explicit analytic form [10] for sim- 
plicity, but it is important to say that Gg^ ilIa*) 
contain some LECs that will represent the main 
source of uncertainty for L\q . 

3. Determination of Effective Couplings 

We will use the recent ALEPH data on 
hadronic r decays [11], that provide the most pre- 
cise measurement of the V—A spectral function. 

The relations ([2]) and ([3]) are exactly satisfied 
only at Sq— >oo, but we are forced to take finite 
values of sq neglecting in this way the rest of the 
integrall. From the SQ-sensitivity of the effective 
parameters one can assess the size of this theoret- 
ical error (quark- hadron duality violation -DV-). 

In Fig. [l] we plot the value of L\q obtained for 
different values of sq , with the one-sigma experi- 
mental error band, and we can see a quite stable 
result at so > 2 GeV^ (sohd lines). The weight 
function 1/s decreases the impact of the high- 
energy region, minimising the DV; the resulting 
integral appears then to be much better behaved 
than the sum rules with s" {n > 0) weights. 



^Equivalently, we are assuming that the OPE is a good 
approximation for n(s) at any |s|=so, what is not expected 
to happen near the real axis and that produces the DV. 
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Figure 1. L^q (sq) and Gff (sq) from different sum 
rules. For clarity, we do not include the error 
bands associated with the modified weights. 



There are some possible strategies to estimate 
the value of L°q and his error. One is to give 
the predictions fixing sq at the so-called "dual- 
ity points" , two points where the first and second 
Weinberg sum rules (WSR) [16] happen to be sat- 
isfied. In this way we get L°q — —(6.50 ± 0.13) • 
10^'^, where the uncertainty covers the values ob- 
tained at the two "duality points" . 

If we assume that the integral ([2]) oscillates 
around his asymptotic value with decreasing os- 
cillations and we perform an average between the 
maxima and minima of the oscillations we get 
L^ff = -(6.5 ±0.2) • 10-3. 

Another way of estimating the DV uses appro- 
priate oscillating functions defined in [17] which 
mimic the real quark-hadron oscillations above 
the data. These functions are defined such that 
they match the data at ~3 GeV^, go to zero with 
decreasing oscillations and satisfy the two WSRs. 
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We find in this way L"^ = -(6.50 ± 0.12) • lO'^, 
where the error spans the range generated by the 
different functions used. 

Finally we can take advantage of the WSRs 
to construct modified sum rules with weight fac- 
tors w(s) proportional to (1 — s/sq), in order to 
suppress numerically the role of the suspect re- 
gion around s ~ so [18]. Fig. [T] shows the re- 
sults obtained with 'Wi{s) = (1 — s/sq) /s (dashed 
line) and 'W2{s) = {l — s/so) /s (dot-dashed line). 
These weights give rise to very stable results 
over a quite wide range of sq values. One gets 
Ll^ = -(6.51 ± 0.06) • IQ-^ using wi(s) and 
Ll^ = -(6.45 ± 0.06) • 10-3 using W2{s). 

Taking into account all the previous discussion, 
we quote as our final conservative result: 

L^ff = -(6.48 ±0.06) • 10"^ (5) 

We have made a completely analogous analysis to 
determine Cfy. The results are shown in Fig. [TJ 
The solid lines, obtained from Eq. ([3]), are much 
more stable than the corresponding results for 
LIq, due to the 1/s^ factor in the integrand. The 
dashed and dot-dashed lines have been obtained 
with the modified weights W3(s) = ^ ^1 — ^^-^ and 

W4(s) =^ (l-^)' (l + 2^). The agreement 
among the different estimates is quite remarkable, 
and our final conservative result is 

= (8.18 ±0.14) • 10^3 GeV^^ (6) 

4. Determination of L\q and C^^ 

The xPT coupling U[q{^) can be obtained from 
L"^, using the relation At 0{p'^) the deter- 
mination is straightforward and one gets 

L[f){p = Mp) = -(5.22 ± 0.06) • 10"^ . (7) 

At order p^, the numerical relation is more in- 
volved because it gets small corrections from 
other LECs. It is useful to classify the 0{p^) con- 
tributions through their ordering within the 1 /Nc 
expansion. The tree-level term Gqj^{^) contains 
the only 0{p^) correction in the large- A^c limit, 
4m^(Cgi — Ci2 — Cgg), that is numerically small 
because of the suppression and can be esti- 
mated with a moderate accuracy [5,8,15,19,20]. 



At NLO GQ^(/i) contributes with a term of 
the form ^1^(0^2 — C13 — Csi)- absence 
of information about these LECs we will adopt 
the conservative range [Cgj- C^^ — Cg^l < [Cg^- 
Cl2—CgQ\/3, which generates the uncertainty that 
will dominate our final error on L^q. Also at 
this order in l/Nc there is the one-loop correc- 
tion G\j^(^) that is proportional to Lg, which is 
better known [21]. Calculating the l/N^ sup- 
pressed two- loop function G^^{ii) and taking all 
these contributions into account we finally get the 
wanted 0{p^) result: 

Lw{Mp) = -(4.06 ±0.04ic«±0.39LECs)- 10-3 
= -(4.06 ± 0.39) • 10-^ (8) 

where the error has been split into its two main 
components. Repeating the same process with 
Cgy (where the only LEG involved is Lg) we get 

C^riMp) = (4.89 ± 0.19) • 10"^ GeV"^ . (9) 

5. Summary 

Using general properties of QCD and the mea- 
sured V — A spectral function [11] we have de- 
termined the chiral LECs L^^iMp) and C^^{Mp) 
rather accurately, with a careful analysis of the 
theoretical uncertainties. 

There are other determinations of Liq from 
r data in the literature. Our result for L\q 
agrees with [12,13], but our estimation includes 
a more careful assessment of the theoretical er- 
rors. The 3.2 cr discrepancy between the estima- 
tion of ref. [14] and ours is caused by an underes- 
timation of the systematic error associates with 
the duality-point approach used in that reference. 
In [13] also C|y is determined with a good agree- 
ment with our result again. The extraction of 
L\^{p?) from has only been done previously 
in ref. [12], at 0{p'^). 

Our determinations of L^Q(Afp) and Cl'j{Mp) 
agree within errors with the large-A^p estimates 
based on lowest-meson dominance [4, 6, 15, 22] 
Lio « -3/2/(8M2) « -5.4 • 10-3 and C^i « 
7/2/(32M^) K. 5.3 • 10-3 GeV-2 and with the 
result of ref. [7] for Cgy, based on Fade Approxi- 
mants. These predictions, however, are unable to 
fix the scale dependence which is of higher-order 
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in I/Nq- More recently, the resonance chiral the- 
ory Lagrangian [6, 23] has been used to analyse 
the correlator n(s) at NLO order in the 1/Nc 
expansion. Matching the effective field theory de- 
scription with the short-distance QCD behaviour, 
the two LECs are determined, keeping full con- 
trol of their /i dependence. The theoretically pre- 
dicted values LlaiMp) = -(4.4 ± 0.9) • 10"^ and 
C^^iMp) = (3.6±1.3)-10-3 GeV-2 [9] are in per- 
fect agreement with our determinations, although 
less precise. A recent lattice estimate [24] finds 
LIq{Mp) = -(5.2±0.5)-10"3 at order p'', in good 
agreement with our result ([7|). 

Using the results of ref. [25], the SU(2) xPT 
LEG I5 can be extracted from L\q{^). We find 
I5 = 13.30 ± 0.11 at 0{p^) and I5 = 12.24 ± 0.21 
at 0{p^). 

Recent analyses of the decay 7r+ l^v^ at 
0{p^) have provided accurate values for the com- 
binations L9-I-L10 [20] and I^—Iq [26], that can 
be combined with our results to get Lg{Mp) — 
(5.5±0.4) -10-3 and le = 15.22±0.39 to order 
that are in perfect agreement with refs. [21,27]. 
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